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We show that the celebrated Painleve equations for the Ising correlation functions 
follow in a simple way from the Ward Identities associated with local Integrals of Motion 
of the doubled Ising field theory. We use these Ward Identities to derive the equations 
determining the matrix elements of the product a{x)a{x') between any particle states. The 
result is then applied in evaluating the leading mass corrections in the Ising field theory 
perturbed by an external magnetic field. 
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1. Introduction 

It is well known since the pioneering work of Barouch, Tracy, Wu and McCoy that 
the spin-spin correlation function of the Ising model at zero external field is expressed in 
terms of suitable solution of certain ordinary differential equation, the so-called Painleve III 
equation. Later this result was rederived and generalized via different approaches 
While Sato, Miwa and Jimbo used the theory of monodromy preserving deformations 
of ordinary differential equations, the derivation of Babelon and Bernard was based on 
exact form-factors 0. One of the goals of this paper is to present yet another derivation 
of this classic result, which we believe is somewhat simpler and more straightforward from 
the field-theoretical point of view. It is based on the well known fact that the Ising field 
theory, being a theory of free fermions, possesses an infinite set of local integrals of motion. 
The differential equation for the spin-spin correlation function follows directly from the 
associated Ward identities. The simplest way of derivation along these lines utilizes the 
local integrals of motion of the "doubled" Ising field theory, which consists of two identical 
copies of the Ising field theory, with no interaction between the copies. Various advantages 
of such "doubling" were recognized before, see e.g. 0,0] ■ The doubled Ising field theory 
is equivalent to free Dirac fermion theory, and it is known to have SL{2) Kac-Moody 
symmetry generated by local charges . The derivation based on the corresponding Ward 
identities is presented in Section 4 below. 

The local integrals of motion act on particle states in a simple way, and therefore the 
approach based on the Ward identities applies as well to all matrix elements of the form 

{AiP[)...A{/3'M)\ a(x)a(x') | A(/3i) . . . ) , (1.1) 

where A{Pi) stand for the Ising free fermions with rapidities Pi, and a{x) denotes the spin 
field. We derive the equations determining these matrix elements in Section 4. 

Our particular interest in the matrix elements ( |1 . 1| ) was motivated by the applications 
to the Ising field theory with nonzero external field h. Adding the term h J a{x) (Px to the 
free-fermion action gives rise to an interesting interacting field theory 0] which describes 
the scaling domain of the Ising model in a magnetic field. Generally, this field theory is 
not integrable, hence no exact solution is available. However, at small field, h « m}^/^, 
perturbation theory in h can produce useful results. As usual, such perturbative calcula- 
tions involve matrix elements of the type ( |1 . IJ ) . For instance, perturbative calculation of 
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the vacuum energy requires the correlation functions, and numerical solution of the corre- 
sponding exact differential equation yields the coefficients of the /i^ expansion with high 
precision 11^,0. The matrix elements ( p. . 1|) involving particles are used in the evaluation 
of the perturbative corrections to particle masses and scattering amplitudes. In particular, 
the ~ h'^ correction to the mass is obtained by evaluating the self-energy part, i.e. 

5m^ = -h^ J { MP) I cr{x)a{0) \ A{/3) )irred d^x , (1.2) 

where (. . .)irred denotes the connected one-particle irreducible matrix element, and the 
integration (which makes the result independent of (3) is taken over the Euclidean plane 
IR^. The integral here can be computed numerically, using exact differential equation de- 



termining the matrix element involved. The results were previously quoted in [|T2[ without 
derivation. Here, in Section 5 below, we fill this gap by presenting some details of this 
calculation. 



2. The Ising Field Theory 

In this section we describe the structure of the Ising field theory in some detail. This 
is done mostly in order to fix our notations and conventions; the reader familiar with the 
subject can safely skip this Section (except maybe for the important Eqs. ( |2.16| ), ( p.l7| )). 

As is well known (see e.g. 0), the Ising Field Theory with zero magnetic field is a 
free Majorana fermion theory, described by the standard euclidean action 

Here we have assumed that the theory is defined on an infinite plane IR^, whose points x 
are labelled by cartesian coordinates (x, y) = (x(x), y(x)), and d^x = dxdy; In what follows 
we typically interpret y as the euclidean time direction. Complex coordinates are defined 
as z(x) = X -f zy, z{x) = x — iy, and the derivatives 9, d in ( p.l|) stand for = ^{d^ — idy) 
and 9z = |(9x + idy), respectively. The chiral components '(/'(x), '^(x) of the Majorana 
fermi field obey the linear field equations 

di^='-^^, d^ = --^i^, (2.2) 

and their normalization in the action (^7l|) corresponds to the following short-distance limit 
of the two-point operator products 

(z — t!) 'i/'(x)'i/'(x') — > 1 , (z — z') ^{x)^{x') — > 1 , as |x — x'| ^ . (2-3) 
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Typical observables of the Ising field theory are described in terms of the order and 
disorder fields, a{x) and fx{x) respectively []T^. These fields are non-local with respect to 
the fermi fields; the operator products 

tp{x)a{xQ) , tp{x)a{xo) , tp{x)n{xo) , tjj{x)n{xo) , (2.4) 

acquire a minus sign when the point x is taken around xq. This property does not define 
the fields a and n uniquely, but a precise definition can be given in terms of the "radial 
quantization" of the theory ( |2.1D . Since the products ( |2.4D as functions of x obey the field 
equations (p.2|) , it is useful to introduce a complete set of solutions, 

Unix) J ^2/ ^ \-ie'^^+-^^'l^^^{mr) 



(2.5) 



(here n E r, 9 are polar coordinates around xq, i.e. z — zq = re*^, z — zq = re /i/ are 
modified Bessel functions) and write the fields -(/'(x) and ifjix) in the products ( p.4|) as 



oo 



(2-6) 

'(/'(x) = ^ anU-n{x) + anV-nix) . 
n= — oo 

The coefficients and a„ here are understood as operators acting on the space of fields; 
they can be easily shown to satisfy the canonical anticommutation relations 

{an,an'} = Sn-i-n',0 , {<in, <in'} = Sn+n' ,0 , {«n,«n'} = 0. (2.7) 

The fields a and n are defined as the "primary" fields with respect to the above alge- 
bra (|2.7|). In particular, they satisfy 



ttn a{x) = , ttn a{x) = , n{x) = , tt^ fJ.{x) = (2.8) 
for n > 0, as well as 

ao a(x) = ^ fi{x) , ao n{x) = ^ a{x) , 

i ^ (2.9) 

ao o-{x) = -j= ii{x) , ao /u(x) = --j= a{x) , 
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where uj = e*^/'^ and uj = e ^'^Z^. These equations do not determine the normahzation of 
the fields cr and /U; we fix it through the short-distance hmit of their operator products 



\x — x'P^^ a{x)a{x') 1 



X 



ctS X X 



. (2.10) 



Applying the operators and dn with negative n to the "primary" fields a and fi 
creates an infinite tower of "descendent" order and disorder fields. It will become important 
below that the lowest of these "descendents" are expressed in terms of the coordinate 
derivatives of the "primary" order and disorder fields themselves. Indeed, for any local 
field O{xo), its derivatives dzoO{xo) and dzoO{xo) can be expressed as 



dz^O{xo) = 7^ / [Tix) dz{x) + e(x) dz{x)] O{xo) , 
d-z,0{xo) = (I \f{x) dz{x) + e(x) dz{x)] 0{xq) 



(2.11) 



where the contour integrals in the x-plane are taken around the point xq in counter- 
clockwise direction and 
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(2.12) 



are related to the components of the energy-momentum tensor and normalized according 
to the usual conformal field theory convention, T = —2nTzz, O = 2nTzz and T = —2nTzz- 
Of course, the contour integrals ( p. 11 ) are nothing else but the commutators with the 
energy-momentum operators P and P, 



dz,OM = [P, O{xo)] , d-z,0{xo) = i [P, O(xo) 



^2.13) 



where 



P = -;^ / [T(x) -t- e(x)] dx 
27r l_r^ 



P--^/ [f (x) -f e(x)] c^x , (2.14) 



with the integration taken along an "equal-time" slice y = constant. If the local field O is 
either a or /i or one of their descendents, it is straightforward to express these commutators 
in terms of the mode operators a^, dn in (|2.6|) , 



^ CO 

n=0 



(2n -f l)a-n-ian + 



n=0 



(2n -f l)a-n-ian + 



2n+ 1 



m 
2n+l 



O'-nO-n+l 



O{xo) , 
O{xo) . 



(2.15) 



Using the defining equations for a and fi, (|2.8| ) and (|2.9|) , tlie above equations ( p.l5|) yield 
the lowest descendents 



a-i a{x) = 4d^{x) , a_i ^{x) = 4da{x) 



72 "-ipv-;-^^ 

a-ia{z,z) = —=4d^{x) , a_i = ^p49a(a;) 
V 2 v2 



(2.16) 



and 



a_2 = , a_2 fi{x) = 



(2.17) 



Of course, the fields a{x) and ^{x) are not mutually local; again, the product a{x)fx{x') 
picks up a minus sign when x is brought around x'. In fact, this product admits operator 
product expansion which involves the Majorana fields '0, t/j and local composite fields built 
from them. Thus, the leading short-distance terms are 



a{xMx') ~ ^ L (z - z'y/' ^(x') + CD (z - zy/2 Vi(x')l (2.18) 

as I X X I ^ . 

The Majorana theory ( |2.1|) describes the critical behaviour of the Ising model in the 
scaling regime near its phase transition point T^. The theory applies to both phases, just 
below and just above the Tc, depending on the choice of the sign of the mass parameter m 
in (|2.1| ). Our definition in ( p.9| ) corresponds to the identification of the case m > with 
the ordered phase T < while the case m < is identified with the disordered phase 
T > Tc. Of course, the theories at the oposite signs of m are related to each other through 
the well-known Kramers- Wannier duality ^ (V', — V')? cr ^ A*? ^^'^ there is no need 

to consider them separately. For this reason in what follows we assume that m is positive, 
i.e. we interpret our results in terms of the ordered phase. In this case it is the order field 
cr which develops nonzero vacuum expectation value [ p^ , 

( a{x) ) = a = m^/s s , 8 = 21/12^-1/8^3/2 ^2.19) 

{A = 1.28243 ... is Glaisher's constant), while the disorder field /i does not, ( /i(x) ) = 0. 

Also, ( pT|) is a free field theory, and of course its Hilbert space of states H associated 
with an "equal-time" slice y = constant is the Fock space of free fermions. It is generated 
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by canonical fermionic creation and annihilation operators, sl^{(3) and a(/3), subject to the 
canonical anticommutation relations 



{at(/3),a(/3')} = 27r5(/3-/3') , {at(/3), at(/3')} = , {a(/3), a(/3')} = . (2.20) 



Here /3 is the usual rapidity variable describing the energy-momentum state of the particle, 
= (m cosh /3, m sinh /3) . The canonical operators are the Fourier modes of the 
Heisenberg field operators if^ix) and i{j{x), i.e. 




— / 7TI 

V'(^) = A/y 



°° df3 

-oo 

°° dp 



^l^'j gy cosh f3 — ix m sinh /3 



ia(/3) e 



-y m cosh /3+ix m sinh /3 



-13/2 



-oo v27r 
Below we use the notation 



y m cosh f3 — ix m sinh /3 _j_ ^ a(/3) C"'^ cosh /3+ix m sinh /3 

(2.21) 



A{(3i)...A{(3N)) = Si\(3i)...a.\(3N)\0) 



(2.22) 



for an A^-particle state. 



3. Conserved Charges of the Doubled Ising Field Theory 

It was observed in the past that many aspects of the Ising field theory simplify if 
one considers a system of two identical copies of the Ising model, with no interaction 
between the copies - the "doubled" Ising field theory [^J^]. Here we also use this trick, 
introducing two copies henceforth referred to as "copy a" and "copy b" . We will use the 
subscript a or 6 to distinguish between the fields belonging to the corresponding copy; thus, 
ifjaix) , ipa{x) , aa{x) , /ia(x) wiU staud for the fermionic and the order- disorder fields from 
the copy a, and ^tpbix) , ipb{x) , ab{x) , ^b(x) will denote these fields from the copy b. The 
correlation functions of the doubled theory factorize in terms of the correlation functions 
of the individual copies, in an obvious manner. 

The above straightforward definition has an inconvenient feature that the fermi fields 
belonging to different copies commute. This can be fixed by introducing, in addition to the 
two Ising copies, two auxiliary variables, the so-called "Klein factors", rja and rj^. These 
are assumed to commute with all the observables from the original Ising copies, and to 
satisfy the following defining relations 

= 1 : ^6=1. VaVb = -VbVa , (3.1) 
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and 

iVa) = , ivb) = . (3.2) 
We then modify the definitions of the fields of the doubled Ising field theory as follows, 

(3.3) 

and 

aa(x) aa{x) , ab{x) -> ab{x) , 

(3.4) 

fJ'aix) ^ rja fJ.a{x) , jlb{x) ^ Tjb Hb{x) ■ 

This modification has no effect on any of the relations in Section 2, and hence it 
does not bring any change to the correlation functions involving only the fields from a 
single copy, either a or h. On the other hand, the Klein factors give rise to additional sign 
factors when factorizing generic correlation functions with entries from both copies present. 
For example, {aa{x)ab{x) aa{x')ab{x')) = (a"(x) cr(x'))^, but {jia{x)jib{x) fXa{x')fXbix')) = 
— {ji{x) fx{x'))'^, where the expectation values in the right-hand sides are that of a single 
Ising field theory. The desirable effect of this modification is that fermi fields from different 
copies now anticommute inside the correlation functions, i{jaix)iljb{x') = — 'i/'b(x')'i/'a(ic), etc. 
Of course, with this definition the doubled Ising field theory is identical to the theory of a 
free Dirac fermi field 

^f(x) = {i^aix) + iifjbix), iJaix) + iif^bix)) . (3.5) 



The space of states of the doubled Ising field theory involves two species of free fermi 
particles, A and B, associated with the copies a and b, respectively. The multiparticle 
states are generated by the corresponding creation and annihilation operators aL^{(3) , a(/9) 
and bt(/3),b(/3): 

\A((3i) . . . A(/3n)B{(3n+i) . . . B{[3m)) = at(/3i) . . . {(3N)h\(3N+i) • • • h^Ml^) • (3.6) 

The operators a"'"(/3) , a(/3) and b"'"(/3) , b(/3) are related to the fields '(pa{x),^a{x) and 
'0b(x), 'i/'f,(x), respectively, through equations identical to Eqs. ( p.21|) . Of course, each pair 
a"'"(/3) , a(/?) and b"'"(/3) , b(/3) obeys the canonical anticommutators ( |2.2C1| ), and a"'"(/3) , a(/3) 
anticommute with b^(/3) , b(/3). 
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Being a free fermion theory, the doubled Ising field theory exhibits an infinite set of 
local Integrals of Motion (IM) . We have no need to describe the whole set here. Certain 
subset of these IM forms the commutator algebra SL{2) of the level zero; its detailed 
description is given in Ref. . Here we display only the fundamental generating elements 
of this subset, since these basic IM are directly used in the calculations below. 

The simplest IM is the U{1) charge associated with phase rotations of the Dirac 
field (pD, 

Zo = — / ['Pa'tpb - ''Pa'tpb] dyi . (3.7) 

J — oo 

Also, the energy-momenta (P, P) associated with each copy are conserved separately, 
hence their linear combinations 

= Pa — Pfe , X_i = Pa — Pb , (3.8) 

are conserved as well. Less trivially, the following integrals are conserved. 



2 



7TI — 

I/Jadl/Jb + i — i^ai>b 



-oo 
oo 

— — — f 1 1 — 



' _ 



(3.9) 



as one can check by a simple direct computation. Equally straightforward calculation 
yields the basic commutators 



[Zo,Xi] = 2zYi , [Zo,Yi] = -2zXi , 

[Zo, X_i] = 2z Y_i , [Zo, Y_i] = -2z X_i 



[Xi,Y_i]=2z(^yj Zo, [X_i,Yi] = 2z(^-j Zq 
[Xi,X_i]=0, [Yi,Y_i]=0, 



(3.10) 



which show that the operators Zo,X±i, Y-ti generate the affine Lie algebra SL{2) of the 
level zero 0. Further commutators of these basic generators give rise to other elements of 



^ The standard Chevalley generators {E±, F±, H±} for the SL{2) algebra are the hnear com- 
binations 

E± = (Xi ±iYi)/m , F± = (X_i^iY_i)/m , H± = ±Zo . 

One can check validity of the Serre relations [E±, [E±, [E±,E^]]] = and [F±, [F±, [F±,F^]]] = 
by yet another straightforward computation. 
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SL{2), which aU have the form of integrals of local densities, quadratic in terms of the 
fermi fields and their derivatives ISll. In our calculations below two such elements, 



2i 



[Xi,Yi 



2iZ. 



[X_i,Y_ 



(3.11) 



will be particularly useful. Explicitly, 

1 



Z2 = 

Z_2 = 



27r 
1 

2^ 



-oo 
oo 



m 



m 



dx , 



(3.12) 



In the next Section we study the Ward identities associated with this SL(2) symmetry 
of the doubled Ising field theory. Derivations of the Ward identities involve commutators 
of the symmetry generators with local fields. Specifically, we will be interested in the 
commutators of the generators with the products of the order or the disorder fields from 
different copies, i.e. with the fields aa{x)ab{x), aa{x)fj,b{x), fj,a{x)ab{x) and /ia(x)/if,(x). 
Note that the densities of the IM ( p.7|) - (|3?9D , ( |3.12|) are local with respect to such product 
fields, and therefore these commutators are again local fields. Of course, the commutators 
with X±i can be written down immediately from the definitions (|3.8|), 



(3.13) 



[Xi, Oa{x)Obix)] = I dOa{x)Ob{x) - i Oa{x)dOb{x) , 
[X_i, Oa{x)Ob{x)] = -Z dOa{x)Ob{x) + I Oa{x)dOb{x) 



Derivation of the other commutators are less elementary, but yet straightforward. One can 
use the decomposition (|2.6|) for each of the fields and i/^b to express any such commutator 
in terms of the corresponding radial mode operators a„, a„ and bmbn- For instance, 

oo 

[Zo, OaOb] = -i {aobo + dobo) OaOb - i ^ {a-nbn - b-nan + d-nK - b-ndn) OaOb . 

(3.14) 

When the primary order or disorder fields are concerned, Eqs. (|2.8D , ( p.9|) , and ( |2.16|) , 
( ^.17| ) apply. It is important to keep track of the Klein factors Ty^, r]b in these calculations 
to get the signs right. For the commutators involving Zq and Yi, the results are 



[Zo,aa{x)ab{x)] = 



[Zo, Ha{x)nb{x)] = 



[Zo, aa{x)fXb{x)] = -i ^a{x)ab{x) , [Zq, fXa{x)ab{x)] = i aa{x)^b{x) , 



(3.15) 
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and 

[Yi, CTa(x)cTb(x)] = -dHa{x)Hb{x) + Ha{x)dHb{x) , 

[Yi, ^a{x)fXb{x)] = -daa{x)ab{x) + aa{x)dab{x) , 
[Yi,(Ja(x)/ib(x)] =id^ia{x)crb{x) -i^ia{x)dab{x) , 
\Yi, lia{x)ab{x)] = -idaa{x)fXb{x) + i aa{x)dfXb{x) . 
Let us also present two useful commutators involving the operator Z2, 

[Z2, aaix)ab{x)] = 2dnaix)dnb{x) - d"^ Ha{x)iJ.b{x) - iJa{x)d'^ Hb{x) , 
[Z2, iiaix)iib{x)] = 2daa{x)dab{x) - 9^cra(x)a{,(x) - aa{x)d'^ab{x) . 



(3.16) 



(3.17) 



The operators Y_i and Z_2 satisfy relations similar to (|3.16| ) and (|3.17| ) with d replaced 
by d and i replaced by — i 

In addition to the above commutators, the calculations in Section 4 below will involve 
the action of the generators on the multiparticle states | A(/3i) . . . A{[3n)B{I3n+i) ■ ■ ■ B{(3m)) 
of the doubled Ising model. The equations 



^0 = 



.t 



(/3)b(/3)-bt(/3)a(/3) 



dp 
J 2^ 



(3.18) 



and 



Yi = ^ 
2 



-00 
00 



at(/3)b(/3)+bt(/3)a(/3) 
t(/3)b(/3)+bt(/3)a(/3) 



2^ 



(3.19) 



follow directly from ( |3.7| ) , ( |3.9D and the decompositions ( ^^.21] ) of the fermi fields i/^a and ■i/'6 
in terms of the particle creation and annihilation operators a, a^ and b, b^. 



4. Correlation Functions and Particle Matrix Elements 

In this section we use the Ward identities associated with the SL(2) symmetry ( |3.1C1| ) 
of the doubled model to derive the differential equation determining the matrix elements 
of the type (|1.1D, involving order and/or disorder fields, between any particle states. As 



the warm-up exercise, let us first rederive the celebrated equations of Ref. |jl|] for the 
vacuum-vacuum matrix elements. 
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4.1. The Correlation Functions 

Consider first the following simple identity 



(0| [ Z2 , aa{x)ai,{x) ^ia{0)^b{0) ] |0) = . (4.1) 
Evaluating the commutator according to Eq. (|3.17|) , one derives the equation 



dGdG-Gd^G-dGdG + Gd^G^O, (4.2) 

where 

G{x) = { a{x)a{0) ) , G{x) = { fi{x)^{0) ) , (4.3) 

and d = dz . Note the signs of the last two terms, which are due to the presence of the 
Klein factors in the definitions ( |3.4| ) of the fields fXh. The equation similar to ([4.2|), 
with d replaced by 9 = 9z , is derived when one takes Z_2 instead of Z2 in ( ^J]) . 

Additional equations are derived in similar way. Again, since the vacuum state |0) is 
annihilated by all the generators (p78|), (|3.9|), the following equations hold 



(0| [X_i,a„(x)a6(x)] [Yi, ^„(0)^b(0)] - [Yi, aa{x)ab{x)] [X_i, ^„(0)^;,(0)] |0) 

2 

= z — (0| [Zo, aa{x)ab{x)] ^ia{0)^ib{0) |0) , (4.4a) 
(0| [X_i,aa{x)^b{x)] [Yi,/x„(0)(Tb(0)] - [Yi,a,(x)//f,(a:)] [X_i, //„(0)ab(0)] |0) 

2 

TTl 

= i — (0| [Zo, aa{x)^b{x)] ^ia{0)a,{0) |0) , (4.46) 

where Zq in the right hand sides comes from the commutator [X_i, Yi] (see Eq. (|3.10|) ). 
Now using Eqs. ( |3.13| ) and ( |3.16| ) one derives 

GddG-dGdG + GddG-dGdG = , (4.5a) 
GddG-dGdG + GddG-dGdG= (^^Y GG . (4.56) 

The Eqs. (^]^) and ( [4.5a, fcj ) are known as the quadratic form |]T^JT^] of the famous differ- 
ential equations of Wu, McCoy, Tracy and Barouch [Q. When G{x) and G{x) are written 
in terms of auxiliary functions f{x) and x(^) ^-s 

m-^/^ G = e^/^ cosh(<^/2) , m"^/* G = e^/^ sinh((/?/2) , (4.6) 
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the Eqs. ( [4 .21) , ( [4.5a, b ) take the sinh-Gordon form, 



m 



ddip = — sinh(2(^) , + (Qipf = , 

a% = - fl - cosh(2v.)l - ^^X+(av^)' = 0. ^'-^^ 



The correlation functions (|4.3|) are related to the rotationally-invariant solution (p = (p{r), 
X = x(^)7 where r = |a;| is the distance between x and the origin; thus, Eqs. ( [4.7| ) reduce 
to ordinary differential equations with respect to r. The solution relevant to the problem 
can be singled out by specifying its r — asymptotic behaviour, 

771 7* 1 

(p{r) ~ - In In (-0) + O (r^ Q"^) , x{r) ~ - In (4m r) + In (-0) + O (r^) ; (4.8) 

here and below in the Appendix we use the notation 

= In (—mr) , (4.9) 



where 7 is the Euler's constant. This solution has the property that <f{r) decays at r ^ 00, 
while x(r) approaches a constant 41ns. All these properties are elaborated in Ref. 
where more details can be found. We present some useful expansions in the Appendix. 

4-2. One- Particle Matrix Elements 

It is straightforward to extend this technique to the matrix elements of the 
form ( p. . 1|) involving the particle states. First, consider the one-particle matrix elements 
(0 I a{x)ii{x') I A{(3) ). It is useful to introduce the following notations 



(0 I a{x)fi{x') I A{P)) = E{x + x'- p) F{x - x'; /3) 
(0 I ^(x)a(x') I A{(3)) = E{x + x'; f3) F{x - x'; (3) 



(4.10) 



where 

E{X + X'- /3) = e-(y+y') ^ cosh/?+z(x+x')f sinh/3 

are the "center of mass" plane waves. The fact that F and F depend on the separation x—x' 
only is a simple consequence of the kinematics. Likewise, from the rotational symmetry 
of the problem, both functions depend only on r and the combination d = 6 — i(3, where 
(r, 9) are the polar coordinates associated with the separation x — x', i.e. 

z-z' = re^^ z-z' = re-^^ (4.12) 
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And it follows from what was said in Section 2 about the local properties of the order 
versus the disorder field that F{x — x'; (3) and F{x — x'; (3) are double- valued functions of 
the coordinates, changing the sign when x is brought around x' . With this understood, 
the above definition implies F{x; P) = ±F{—x; (3), where the sign depends on the way of 
continuation x —x. 

The differential equations for the functions F and F can be derived from the identities 



( I Yi aa{x)ai,{x) a„(x')/Ub(x') | A{(3) ) = , 
( I Yi aa{x)^b{x) a„(x>b(x') | A{(3) ) = . 



(4.13) 



One uses the commutation relations (|3.16 ) to move the operator Yi to the right, and 



then applies (|3.19|) to evaluate the action of this operator on the particle state. The 
resulting equations are linear in F and F, and involve the correlation functions ( [4 .31) as 
the coefficients. They are brought to a nice form by introducing the functions and 
related to F and F as 

F + iF = ivm^/^e^/^^+ , 



(4.14) 



where again uj = e*'^/^, uj = e '^^/^^ and x is related to the two-point correlation functions 
as in ( [4.61 ). With these notations, the Ward identities ( [1.13D yield 



+ 



m 



(4.15) 



+ ? 



where d = dz and the function (p is the same as in ( [4.6[ ). Replacing Yi in ( |4.13D with Y_i 
one derives similar equations involving the derivative d = dz , 



(4.16) 



One recognizes in ( [4.15| ), ([4.16| ) the Lax representation of the sinh-Gordon equation: 
the first of the Eqs. ( [4.7| ) guarantees compatibility of (|4.15D and ( [4.16| ). Thus, the one- 
particle matrix elements ( [4.1C1| ) are related to the solution of the linear problem ( [4.15| ), 
( [1.16|) associated with the sinh-Gordon equation (|4.7| ). 
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Let us briefly describe here some elementary properties of the functions '^±. As was 
mentioned above, they depend on two variables, r and d = 6 — iP, where (r, 9) are the 
polar coordinates deflned in ( ^.121 ). The monodromy properties 



^'+(r, ^? + 7r) = z\I/+(r, ^?) , \&_(r, ^? + 7r) = -z\I/_(r, ^?) 



follow directly from the definitions (|4.10|) , (|4.14|) . Also, it is possible to show that 



(4.17) 



(4.18) 



and that when both r and are real (i.e. the rapidity P is continued to pure imaginary 
values) and take complex conjugate values. The r —>■ behaviour of and 
follows from the short-distance operator product expansion ( p.l8|) . 



exW/2^±(r, i9) ~ V2^{mry/^e^'^/^ . 



(4.19) 



More details on the asymptotic behaviour is presented in the Appendix. Note that in view 
of (|4.17| ), (|4.18| ), \l/+(r, t?) and \E'_(r, t9) admit the Fourier decompositions of the form 



(4.20) 



where \l/2n(^) are real at real r. The coefficients \E'2n(^) can be interpreted in terms of the 
finite-distance version of the operator product expansion ( |2.18| ) , i.e. 



a{x)fi{x') = 



n=0 



+ 



(4.21) 



where the dots stand for the contributions of the multi-fermion operators (i.e. the com- 
posite fields built from three or more fermi fields, defined in such a way that their matrix 
elements between (0 | and | A{P)) vanish). Using the polar coordinates ( [1.1 2|) to represent 
the separation x — x' , we have 



(4.22) 



where 



C2n(r) 
C2n-l{r) 



( \n n2n— 1 

^X{r)l2 ^ , 



(-)-2 



n o2n — 2 



(4.23) 



7rm 



2n-3/4 
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4-3. Multi-Particle Matrix Elements 



The matrix elements involving more then one particle can be determined in a similar 
way. Here we explicitly elaborate the matrix elements 



G{x, x';/3i,..., P2n) = ( I a{x)a{x') \ A{f3i) . . . A{P2n) ) , 
Gix, x'; . . . , (32n) = ( I iJ.ix)iiix') \ . . . A{(32n) ) ■ 



(4.24) 



In order to avoid lengthy expressions, we shall drop the explicit indication of position 
dependence from now on and use instead the notations G{(3i, . . . , (32n) and G'(/3i, . . . , (32n) 
for the matrix elements ( f4.24| ). Similarly, \l/±(/3) will stand for the functions defined in 
( [4.14|) in relation to the one-particle matrix elements (^4.101 ). 

Let us start with the derivation of the two-particle matrix elements G(/3i,/32) and 
G{(3i,(32)- It suffices to consider the identities 



(0 |Zo aa{x)^ib{x) aa{x')iJb{x') I A(/3i)S(/32) ) = , 
( I Yi aa{x)at{x) a„(x')a,(x') | A{/3i)B{l32) ) = 
( I Yi fia{x)fibix) fiaix')fibix') \ A{(3i)B{P2) ) = 



(4.25) 



Again, one uses the commutators ( p.l5| ), ( |3.16[ ), and then Eqs. ( |3.18D , ( |3.19| ) to evaluate 
the action of the operators Zq and Yi on the two-particle states. The relations that come 
out allow one to find explicit expression for the desired two-particle matrix elements in 
terms of the one-particle ones. The result is 



.G+(A,/32) _ E{P^)Ei(32: 



G + G 



e/3i -|_ e/32 



^ G-G ef^^ + eP^ 



e^ivE'_(/3i)vE'+(/32) - e^^vl/+(A)vE'_(/32) 
e^ivp+(/3i)vl/_(/32) - e^^vi/_(/3i)vl>+(/32) 



(4.26) 



where 



(4.27) 



G±{Pi,(32)=G{(3i,P2)±G{(3i,P2) , 

and E{(3) = E{x + x'; (3) are the plane- wave factors ( |4.11|) . 

Once the two-particle matrix elements are found, all the multiparticle ones are ob- 
tained through the recursive equation 

2N-1 

G G{Pi, . . . , P2n) = YI i-y^'G{f3i, . . . , . . . , P2N-i)G{P,, P2n) , (4.28) 
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which follows directly from the Ward identity 



( I Zo a„(x)a,(x) aa{x')<n,{x') \ . . . A{^2N-i)B{P2n) ) = . (4.29) 

Of course the recursive equation ( [4.28| ) is nothing but the statement that the multiparticle 
elements ( [1.24[) are expressed in terms of the two-particle ones, through all possible Wick 
pairings, 

G^-i . . . , = e«-«-G(/3„, , . . . , , (4.30) 

where £"i -"2iv (denotes usual antisymmetric tensor and summation over = 1, . . . , 2N is 
implicit. This form reflects the free-fermion structure of the Ising field theory. The matrix 
elements G{/3i, . . . , (32n) have the same form in terms of the functions G'(/3i, (32). 



5. Corrections to the Masses 

As was already mentioned in the Introduction, the primary motivation of this study 
was the development of more efficient perturbation theory for the Ising field theory in the 
presence of external field. The last theory corresponds to adding the spin term to the free 
action (PTT|), 

A = Aff + h (fx a{x) , (5.1) 



and it appears to be highly nontrivial interacting field theory. Its physics depends on a 
single scaling parameter ry = m/|/i|^/^^ and shows rich behavior in terms of its particle 
spectrum [S]. Exact solution exists only in the limiting cases ?7 = ±oo, where of course the 



theory reduces back to the free fermions (|2.1| ), and at = 0, where it is integrable, but at 
generic values of this parameter the mass spectrum is understood mostly on a qualitative 



level, as described in P jr7| , p^ . In this situation, perturbative expansions around r] = ±oo, 
say, are useful. Here we will display the calculations of the leading perturbative correction, 
~ /i^, to the particle mass in the disordered regime r] — oo, and of the same order 
correction to the "quark mass" in the ordered regime — +oo. 

Both these calculations are based on the formula (|l]^) , but of course their interpreta- 
tions in terms of the physical mass spectrum are rather different. In the disordered regime 
r] ^ ~oo the spectrum contains a single particle, and the calculation of the mass correction 
according to (|1.2| ) yields just that - the correction to its mass Mi, i.e. the coefficient a in 

Mi(r7) = |m|+a JjjL + o((-r7)-i5/2) , 77 -oo , (5.2) 
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where m is the mass parameter of the unperturbed theory (|2.1j ). Note that we put |m| here 
because in our convention ( |2.9| ) the disordered regime corresponds to negative m in ( |2.1| ). 
RecaU that all equations in the previous Sections are written down under the opposite 
assumption m > 0. But thanks to the duality m ^ — m, a ^ ^ we can still use them in 
this computation, just taking in ( |1.2| ) the matrix elements of the disorder fields fi instead 



of a. Therefore we have for the coefficient a in (|5.2 



a^-^lmr^/^y {f3\ fj.ix)ij{0)\ (3)i,,edd^x . (5.3) 

On the contrary, in the ordered regime ry +00 the mass spectrum is rather complex. 
The last term in ( |5.1j ) creates a long-range confining iteraction between the original par- 
ticles of ( |2.1| ) , which become "quarks" . The particle spectrum consists of a tower of their 
meson-like bound states. The meson masses Mi , i = 1, 2, . . . , do not admit straightforward 
perturbative expansions in /i^, instead they expand in fractional powers of h P,p!2l, 



{2sf^ _ {2sf'\- fUzl 57 \ (2.-)^ 
2 775/4 40775/2 1^2800 560 y 7715/4 



M.(77) = 2m,(r7) 1+ ^ ,1,,, - \/^,/, ^ + - tt^ + ■ ■ ■ } ^ (M 



where s relates to the spontaneous magnetization at zero h as in Eq. ( |2.19| ), and —Zi are 
consecutive zeros of the Airy function, Ai{—Zi) = 0. Here 771^ denotes the "quark mass" 
which admits usual perturbative expansion in powers of /i^. 



mq{r]) =m + aq^jj^ + 0{'r]^^/'^) , 77^+00, (5.5) 



77 

where the coefficient is computed according to the Eq. (|1.2| 



Go = 

^ 2 



^ mi5/4 j {(3\ a{x)a{0) \ f3 )irred d^x . (5.6) 



Note that the leading radiative correction in ( |5.5| ) competes with the last term in ( p.4| ). We 
also would like to stress here that, since the "quarks" do not appear in the asymptotic states 
of the theory, there is no definition of the "quark mass" independent of the perturbation 
theory. Therefore, while (|5.3| ) is derivable as usual from the LSZ theory, the perturbative 
expansion (|5.5[) is rather a definition of the quark mass. Nonetheless, the leading correction 



in (|5.5|) shows up in many places besides the expansion (|5.4|) of the meson masses, notably 
in the 77 — > 00 expansion of the "false vacuum" resonance decay rate (see [1T2|). 
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With all this said, we just proceed with evaluating the integrals in (|5.3|) and (|5.6|). 
The matrix elements involved are obtained from the matrix elements in Section 4 through 
usual crossing equations. One writes the full matrix elements as 

I a{x)a{x') I A{^2)) = 2n 5{(3i - (32) G + G{(3i\^2) , (5.7a) 
\^i{x)^i{x')\A{(32)) = 2-k5{(3^-(32)G + G{(3^\(32) , (5.76) 

where the delta-function terms are the disconnected "direct propagation" parts shown 
schematically in Fig. 1, with G = G{x,x') and G = G{x,x') being the two-point func- 
tions (|4.3|) , and the functions G(/9i|/32) and G'(/3i|/?2) (whose arguments x,x' are also su- 
pressed) are related to the two-particle matrix elements G{l3i, (32) and G{(3i, (32), Eq. ( [4.24|) , 
as 

G{Pi\P2) = G{Pi-in,P2) , GiPi\P2)=Gi(3i-tn,(32) . (5.8) 



a(x) 




a(x') 



Fig. 1: Disconnected "direct propagation" part in ( 5.7a| ); the double line represents 
two-point correlation function G{x,x') and the plain line indicates a particle of 
rapidity (3i. 



The Eqs. ( |5.3| ) and ( |5.6| ) involve only the diagonal elements G{(3\(3) and G{(3\(3). 
From ([4.26| ) we find 



G(/3|/3)=vl>+(/3)vl>_(/3) 
G(/3|/3)=vl>+(/3)vl>_(/3) 



G-G 



d 



In 



^'+(/9) 



d(3 V^-(/3) 
^+(/3) 



d(3 V*-(/3) 



(5.9) 



Remaining disconnected and one-particle reducible parts still have to be subtracted before 
one plugs these functions in (|573|), (|5.6|), 



(5.10) 



' A{P) I ^{x)^i{x') I A{P) )i,,ed = G{(3\P) - Sim , 
'Ai(3) I aix)a{x') \ Ai(3) )i„ed = Gi(3\(3) - Si(3\(3) . 
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Since the matrix elements {A{(3) \ n{x) \ A{(3')) vanish, but 



(0 I n{x) I A{(3) ) = ae 



-y m cosh /3+ix m sinh /3 



(5.11) 



the matrix element of fx{x)fx{x') has no one-particle reducible parts, and one just has to 
subtract the disconnected pieces depicted in Fig. 2, 



Sim = {Am ^(^) |0) (0| fi{x') \A{(3)) + (0| ^(x) \Am {A{P)\ ^{x') |0) 
= 2 o"^ cosh [m r sin t?) , 



(5.12) 



where d = — i(3 and (r, 9) are polar coordinates ( ^.12| ) associated with the separation 
X — x' . On the other hand, ( | a{x) \ A{j3) ) = 0, so that the matrix element of a{x)a{x') 
does not have such disconnected parts, but it has one-particle reducible components shown 
in Fie;. 3. 



^(x) 



(a) 



ij=e-ip 



P 



\i(x') 



(b) 



Fig. 2: Disconnected terms in the matrix element ( ^(/3) | /x(x)/x(a;') | A{(3)). Sum 
of these diagrams, Eq. ( |5?T2D , is the subtraction term S{j3\l3) in ( [5?ToD . 




a(x) 



ij = e-i p 



a(x') 



(a) 




(b) 



Fig. 3: One-particle reducible components in {A{f3) \ a{x)a{x') 
add up to the subtraction term S{j3\(3) in ( ^.10| ) - see Eq. ( |5.13[ ). 



A{f3)). These 
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This reducible part is written as the principal-value integral 



Sim - r f 



{Ama{x) \A{f3')) {A{(3')\ a{x') 
+ (0| a{x) \Ai(3)A{(3')) {A{f3)Ai(3')\ a(x') |0) 



(5.13) 



over the rapidity /3' of the intermediate particle in Fig. 3; the matrix elements involved 
here are the cr-field form- factors H, 



; I a(0) I A{(3i)A{(32) ) = ia tanh 
I (j(0) I A{(32) ) =ia coth 



I3i-I32 
2 

P1-P2 



(5.14) 



and in writing ( |5.13| ) we have assumed that y > y'. For our purposes the following 
equivalent form of this term will be more convenient, 



SiPIP) = 2a^ mr cos^ + a' e 



— 2 „ — m r sin 1 



I — 2 ^mr sini} 



-Ko{mr) + 2d^A{r,'d-Tv/2) 



(5.15) 



where 



A{r,^) = I ^ tanh (^^y^ ) osh/3' 



-00 
00 



B{r,^) = I ^ coth (^^^ ] e— 



-IT <^e'd < TV 



0<^e^<2Tv. 



(5.16) 



Here again d = 6 — iP, and {r,9) are the polar coordinates ( fl.l2|) . The above integrals 
define the functions A(r, 1!}), B{r, i}) in the specified domains of ^? only; outside them, they 
are defined by analytic continuation. 

Using (15.101 ) , the integrals in ( |5.3|) and (|5.6D were evaluated numerically, after finding 
numerical solution of the differential equations ( [4.15|) , ([4.16D for the functions in ( |5.9|) . 
Since the integrals do not depend on (3 we set 



13 = 0., 



(5.17) 



so that d = 6, just the euclidean angle between the separation x — x' and the x-axis. In this 
case (and indeed at any pure imaginary (3) the functions and are complex-conjugate 
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to each other. Few words about integration of the differential equations (|4.15|) , (|4.16 ) are 
worth saying. We found it most convenient to use the representation of the solution in 
terms of the Backliind transformation of the angular-symmetric solution f{r), x(r) of the 
sinh-Gordon system (^4.7| ) appearing in the representation (f4.6| ) of the correlation functions. 
The solution of the linear system ( ^4.15| ), ( [4.16| ) can be written in terms of x(r) and two 
aiixiliary functions (j){r,i!}) and p{r,'d), 



(5.18) 



where cj) and p must satisfy the equations 



m 



cosh ((/? 



— e ^ cosh (^if + , 



(5.19) 



and 



5(p - X) = y e'^ sinh (v? - i 



d{p-x) = ^ sinh {ip + , (5.20) 



and hence also the sine-Gordon system 



m 



dd(j) = sin(2^) 



ddp 



m 



1 + cos(2(/)) 



d^p - {d(pf = , 
d^p- {BcPY = . 



(5.21) 



The solution of ( p.21| ) which appears through ( |5.19| ) is rather interesting. The phase (p is 
not a single- valued function of the coordinates. Instead, when written in terms of the polar 
coordinates (r, t?), it is quasiperiodic function of the angle, ^(r, + 27r) = ^(r, d) + 27r, as 
demanded by the monodromy properties of the matrix elements ([4.10| ) stated in Section 4. 
Qualitatively, it can be described as the juxtaposition of two sine-Gordon domain walls 
(i.e. the sine-Gordon soliton solutions, in the euclidean nomenclature) of oposite sign, 
extending along the x-axis in oposite directions; the solution is singular at r = 0, and its 
shape in the "junction" region, r ~ m~^, is shown in Fig. 4. 
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2ji 



3JI/2 



u 

-e- 



jr/2 




Fig. 4: Plot of the phase (f>{r,'d) for fixed r as function of = 0. It satisfies 
the symmetry properties (f>{r,'d) + 0(r, tt — -i?) = tt and 0(r, i?) = —(f>{r, —"&). The 
"domain walls" at ~ and ~ tt are clearly visible in the plots for r > 3m~^. 

The numerical integration was performed in three steps. First, the ordinary differential 
equations, the radial form of ( [4.7|) , were integrated numerically, the initial conditions at 
small r being provided by the short-distance expansions ( [A.2| ), and the large-distance 
asymptotics ( |A.5| ) were used for precision control. Next, the functions (j){r, -d) and p(r, •&) 
were computed by integrating the first-order differential equations ( p.l9| ), ( |5.20D , with the 
initial conditions fixed using the short-distance expansions ( |A.6D . This step provides the 
numerics for G'(OlO) and G(0|0) as functions of (r, ^?); according to ( |5.9|) , 



Gim 

G{m 



,P-x/2 



,P-x/2 



(5.22) 



smh.{(f/2) — cosh((/?/2) 
cosh((/?/2) — sinh((/?/2) d^c, 

After the subtractions in ( p.lOj) are made to eliminate the reducible parts, the integrands in 
( ^73|) and (|5.6|) decay exponentially at large r, rendering the integrals convergent. So, as the 
last step, the integrations over the euclidean coordinates r, 9 were performed numerically, 
yielding 

a = 10.7619899(1) , (5.23) 

and 

ttq = 0.142021619(1) . (5.24) 

These numbers provide more accurate data for the mass corrections previously estimated 
in ra. 
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6. Discussion 



The purpose of this paper was two-fold. Firstly, we have presented a new derivation 
of well-known result - the nonlinear differential equations of Ref. for the correlation 
functions in the Ising field theory with zero magnetic field. This part has a methodical 
value at the best. Nonetheless, we decided to include it here because we believe our 
derivation is somewhat simpler then the traditional ones, and also it has certain potential 
for generalizations. It is based on the Ward identities associated with the special integrals 
of motion of the doubled Ising field theory. Incidentally, the "doubling" , albeit convenient, 
is not the most essential part of our approach. The single Ising field theory has a system 
of local integrals of motion powerful enough to render similar derivation possible; the 
"doubling" trick just makes it shorter. Anyhow, our approach can be adopted to yield 
simple derivations of other results in the Ising field theory. Thus, it was already used 
in in deriving the finite-size form- factors of the spin operator, and recently in 



in a simple derivation of the differential equations of in the Ising field theory on a 
Poincare disk. Among interesting potential applications let us mention the Ising field 
theory at finite temperature T (equivalently, the theory defined on euclidean cylinder, 
with the points (x, y) and (x, y + T~^) identified). One can notice that all arguments of 
Section 4.1 remain valid when the vacuum expectation values (0 | ... | 0) are replaced 
by the thermal averages tr(e~^/^ • • •)' because the total Hamiltonian H of the doubled 
theory commutes with all the generators Z±2, X±i and Y±i in (|4.1|) and ( [4.4a, It 
follows that the Eqs. (jOD and 

remain valid in this situation as well - the result 
previously obtained by different methods in j^^] and 0. In this case the functions cp 
and X depend on two variables (x, y), but the relevant solution of the Eqs. ( [4.7| ) still can be 
fixed uniquely by imposing the Matsubara boundary conditions ^pix^y + T~^) = {p{x,j), 
as well as appropriate asymptotic conditions at |x| — oo and (x, y) — > (0,0). We plan 
to exploit advantages of our technique in relation to this problem elsewhere. Yet another 
potential application is in the Ising field theory with boundaries. On the other hand, 
extension of this approach to other integrable field theories is more problematic. The fact 
that all descendants of the spin fields up to sufficiently high level are expressed through 
the derivatives, as in the Eqs. (p.l6| ), ( p.l7|) , was very essential for the whole scheme to 
work; this fact seems to be rather specific feature of the Ising field theory. 



^ Additional equation proposed in [21], the Eq. (4.42d), apparently does not hold. 
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Secondly, the results of Section 4 were used in Sect. 5 in computing the leading per- 
turbative mass corrections in the Ising field theory with magnetic field ( ^.Ij) , and from this 
point of view the numbers ( |5.23D and ( |5.24| ) constitute the main result of this work. This 
computation is a part of our ongoing project of systematic study of the field theory ( |5.1|) . 
We are planning to further apply the technique developed here to perturbative calculations 
of other relevant quantities, such as scattering amplitudes and various structure functions. 
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Appendix A. 

Here we present the short- and large-distance expansions of the functions f{r), x(r) in 
( [4.6|) (mostly borrowed from Ref. as well as the analogous expansions of the functions 
^±(^5^) in ( ^4.14| ). Here again d = 6 — ip, where {r,9) are the polar coordinates ( ^.12| ). 
These expansions are used in the numerical evaluation of the self-energy parts ( |5.3|) and 
( ^^) , as explained in Section 5. 



A.l 

The functions (^(r), x(r) in ( [4.6| ) obey the radial form of the Eqs. (|4.7|) , i.e. 

+ ^dr^ = 4- smh(2(/p) , IdrX = {dr^pf + ^ (1 - cosh(2^)) . (A.l) 

As is explained in [||, the relevant solution is characterized uniquely by its short-distance 
asymptotic behavior ([4.81). Corrections to this leading asymptotics can be obtained by 



iterations of ( |A.l 



ip{r) = - In (^— j -ln(-O) + {mrf + {m rf h + O {{m rf'^ Q.^) 



(A.2) 



1 ( 777* V I 

X(r) = -ln(4mr) -Mn(-O) + — + {mr)'^h4 + (mr)^ hs + ((mr)^^O^) , 

2 8 
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where O was defined in ( [4.9| ), and the coefficients fn and hn are rational functions of O, 

/4 = -TJT^(2^^-l) (402-20 + 1) , 

^ (A.3) 

/s = -T^^jFT^ (2048 O^ - 4096 O^ + 3648 O^ - 1568 O^ + 136 O^ + 1110- 32) , 
2 ° O 



and 



/i4 = -7Ti^(803-802 + 20 + l) , 



hs = — (2048 O^ - 4096 O^ + 3776 O^ - 1888 O^ + 496 O^ -1110 + 32) 



(A.4) 



These equations extend the corresponding results presented in ||^. 

Alternatively, the same solution can be characterized by its large-distance behavior 



(^(r) = ^Ko(mr) + 0(e-3"^'') , 

ATT) 7^ r "1 

[K^{mr) - Kf{mr)] + Ko{mr)Ki{mr) +0 (e"^""^) 



x(r) = 41ns — 



7T' 



m r 



where Kjy are usual modified Bessel functions and s is defined in ( |2.19| ). 



(A.5) 



A.2 

The short-distance expansions of are best written in terms of the representation 
( ^■18|) . The functions ^(r, t?), p{r,d) in ( |5.18| ) solve the first-order differential equations 
( ^■19|) , ( |5.20|) , with the leading r asymptotic (j){r,'d) t?, p{r,'d) ^ ln(47r2mr) 
(see Eq. ( [1.19| )). Corrections to this asymptotic are obtained directly, by iterating ( |5.19|) , 
( ^■20|) , using the expansions ( |A.2|) for ip{r) and x(r). We present just few first terms: 



(/)(r, ^?) = "i? + (mr)2p2 + (mr)'*p4 + (mr)^ pq + O ((mr)^ O^) , 

p(r, t?) = - ln(47r2m r) + (m r)^ q2 + (m r)^ ^4 + (m r)^ Q'e + O ((m r)^ O*) 

where again O is the logarithm ( f4.9| ). 



(A.6) 



P2 



-^O sin2^? 
23 



P4 = — 0(40- 1) sin 4^? 
2 



P6 



— [3 (10 O^ - 7 O + 2) sin 2^? - O (16 O^ - 6 O + 1) sin 6^] 
2 3 
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(A.7) 



and 



g2-^ + ^(20-l)cos2i9^ 
54 -^[(160 -8 



(A.8) 



Q6 



1 



[-9 (12 _ 10 O + 3) cos 2t? + (96 - 36 + 6 O - l) cos 6^] 



216 32 

Large-distance behavior of these functions foUows from exact expansions 



-iTrn/2 



F-iF 



a 



n=0 



n=0 



i/2 



ni 



'S^(r, d - 7T/2) + i e-^'^^^"^A„(r, ^ - 7r/2) 



(A.9) 



of the matrix elements ( |4.10|) . When -d lays inside the strip < ^Rci? < tt, the functions 
An{r, t?) and Bn{r, 1!}) here are defined as the 71-fold integrals 



An{r.^)= n 

l<j<n >- 

Bn{r,^)= n 

l<j<n 



2tH 



-00 
00 



tanh 



coth 



I3j - id 
2 



-m r cosh /3j 



-m r cosh /S, 



JJ tanh' 

l<j<A;<n 

JJ^ tanh^ 

l<j<fc<n 



2 / 1-^3 



13, - 13k 



2 / A^J 



(A.IO) 

(Ao(r, 7?) = BQ{r,d) = 1) and analytic continuation of these functions to full complex 
7?-plane is achieved using the following relations 



S^(r, -d) = {-irB^{r,d) - 27i(-l)"e-^'^°^^A„_i(r,7?) 



(A.ll) 



and 



(A.12) 



A^(r, 7? + tt) = S^(r, 7?) - 2ne-^^°^''A^_i(r, 7?) , 

S^(r,7? + 7r) = A„(r,7?). 

The last relations ensure right monodromy properties of the matrix elements ( [4.101 ). The 
expansions (|A.9|) are easily obtained in a standard way, by writing down the intermediate- 
state decompositions of the matrix elements ( [1.1 0|) , and using the exact form- factors of ^ . 
In fact, the series representations ( [A.9[) can be taken as the starting point of alternative 
derivation of the equations ( [1.15[) , ([4.16[) i. 



This was done independently by F. Smirnov (private communication). 
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